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We present the results of searches for B decays to charmless final states involving <f>, /o(980), and 
charged or neutral p mesons. The data sample corresponds to 384 xlO 6 BB pairs collected with 
the BABAR detector operating at the PEP-II asymmetric-energy e + e~ collider at SLAC. We find no 
significant signals and determine the following 90% confidence level upper limits on the branching 
fractions, including systematic uncertainties: B(B° — * 4xf>) < 2.0xl0~ 7 , B(B + — ► <pp + ) < 30xl0~ 7 , 
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B(B° -> <)>p ) < 3.3xHT 7 , B[B° -> 0/ o (98O)] x B[/ (980) Tr+Tr"] < 3.8xHT 7 , and 
/o (980) /o (980)] x S[/o(980) tt+tO x B[/ (980) K + K~] < 2.3xHT 7 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



We report the results of searches for the decays B° — ► 
0/ o (98O), /o(980)/ (980), and B± -> 0p± [1] 
using data collected with the BAS4R detector. The 
B° — > cfxfi decay is an OZI suppressed process with an ex- 
pected branching fraction in the range (0.1 to 3) xl0~ 8 
in the Standard Model (SM) [2-4]. The decays B° -> </>p° 
and B + — » 0p + are pure 6 — ► c? loop processes; the ex- 
pected branching fractions for these modes range from 
(2 to 7)xlCT 8 [5-9]. The presence of new physics (NP) 
would give rise to additional amplitudes that could en- 
hance the branching fractions for these decay modes rel- 
ative to the SM predictions [2, 3, 6]. The branching frac- 
tion for B° — > <f>(j) could be enhanced to 1CP 7 [2], and 
the branching fractions for B — > 4>p decays could be en- 
hanced by 20% [8] in the presence of NP. We are not 
aware of branching fraction predictions for B° — ► <fifo 
and B° - /o/ . 

The B decays to (fxp and (f>p are complicated by the 
presence of one amplitude with longitudinal polariza- 
tion and two amplitudes with transverse polarization. 
The fraction of longitudinally polarized events is de- 
noted by Jl- Integrating over the angle between the 
vector meson decay planes, the angular distribution 
(l/r)rf 2 r/dcos6»idcos6» 2 is 



-~ — ^ « — — — — — — - — 

/ L cos 2 9 X cos 2 2 + 1(1 - f L ) sin 2 ft sin 2 6 2 



, (1) 



where the indices 1, 2 label the two vector mesons in the 
final state, and the helicity angles 6\ : 2 are the angles 
between the direction opposite to that of the B° (B + ) 
and the K + or ir + (tt°) momentum in the <p or P° (p + ) 
rest frame. We define the angles 6\ t 2 for fo mesons in 
an analogous way. The expected values of /l range from 
0.6 to 0.8 [3, 4, 6, 7] for B° 00, 4>p°, and B± 0p±. 
The presence of NP could lead to enhancements of the 
transverse polarization amplitudes [2, 3, 6]. 

The current upper limit on the B° — ► <p<f> branch- 
ing fraction, obtained from a data sample of 82 fb _1 , 
is 1.5 xl0~ 6 [10]. The upper limits on B° -> (f)p° and 
B + — » 0p + , determined using 3.1 fb _1 of data, are 
1.3xl0~ 5 and 1.6xl0~ 5 [11], respectively. Using a data 
sample of 349 fb _1 , BABAR recently reported an upper 
limit of 1.6xl0- 7 for B° -» / / [12]. This last result 
relies on the assumption that the fo — > 7r + 7r~ branching 
fraction is 100%. In this analysis, we make the compli- 
mentary assumption that one fo decays to 7r + 7r~ and the 
other to K + K~ and search for B° — > /o/o in a cleaner 
final state than Ref. [12]. All these limits correspond to 
a confidence level (C.L.) of 90%. 

The results presented here are based on an integrated 
luminosity of 349fb~ 1 , corresponding to (384 ±4) million 



BB pairs. These data were recorded at the Y(4S) res- 
onance with a center-of-mass (CM) energy s/s = 10.58 
GeV. The BABAR detector is described in detail else- 
where [13], and is situated at the interaction region of 
the PEP-II asymmetric energy e + e~ collider located at 
the Stanford Linear Accelerator Center (SLAC). We use 
Monte Carlo (MC) simulated events generated using the 
GEANT4 based [14] BABAR simulation. 

Photons are reconstructed from localized deposits 
of energy greater than 50MeV in the electromagnetic 
calorimeter that are not associated with a charged track. 
We require 7 candidates to have a lateral shower pro- 
file [15] that is consistent with the expectation for pho- 
tons. 7T° candidates are reconstructed from two 7 candi- 
dates with invariant mass 0.10 < m 77 < 0.16 GeV/c 2 . 

We use information from the vertex detector, drift 
chamber and detector of internally reflected Cherenkov 
light to select charged tracks that are consistent with 
kaon or pion signatures in the detector [16]. We re- 
construct <ft (p ) candidates from pairs of oppositely 
charged kaon (pion) candidates with invariant mass 
0.99 < rriKK < 1.05GeV/c 2 (0.55 < < 1.05GcV/c 2 ). 
For p° candidates we require the helicity angles to sat- 
isfy |cos#i| < 0.98 since signal efficiency falls off near 
I cos 6i\ = 1 . Charged p candidates are reconstructed 
from a charged track consistent with the pion signa- 
ture and a ir° candidate. The invariant mass m^o of 
the p + candidate is required to lie between 0.5 and 1.0 
GeV/c 2 . We also require that the helicity angles satisfy 
—0.8 < cos 9i < 0.98 as signal efficiency is asymmetric be- 
cause of the 7T° meson, and falls off near cos 6i = ±1, and 
background peaks near — 1 . We select fo candidates from 
two charged tracks that are both either consistent with 
the kaon or the pion signature in the detector. We ap- 
ply the same selection criteria to fo — > 7r + 7r~ candidates 
as for p° mesons. Similarly we apply the same selection 
criteria to fo — > K + K~ candidates as for (f> mesons as 
the minimum rriKK we can reconstruct in the detector is 
0.99 GeV/c 2 . 

We reconstruct signal B candidates (B rcc ) from com- 
binations of two <j) mesons, one </> and one p or fo, and 
two fo mesons. The /o/o mode is required to have one fo 
decaying into n + n~ , and the other decaying into K + K~ . 
We require the fo in <f>fo to decay into 7r + 7r~. 

We use two kinematic variables, toes and AE, in order 
to isolate the signal: toes = \J ( s /2 + Pi • PB) 2 /Ef — p 2 B 
is the beam-energy substituted mass and AE = E* B — 
y/s/2 is the difference between the B candidate energy 
and the beam energy in the e + e~ CM frame. Here the 
B roc momentum ps and four-momentum of the initial 
state (Ei,pi) are defined in the laboratory frame, and 
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Eg is the -B rec energy in the e + e~ CM frame. The dis- 
tribution of 7Tt.es (AE) peaks at the B mass (near zero) 
for signal events and does not peak for background. We 
require toes > 5.25GeV/c 2 . For the <jxj) final state we 
require \AE\ < 0.15 GcV. To reduce background from 
non-signal B meson decays we apply the more stringent 
cut of -0.07 < AE < 0.15 GcV for all other modes. 

The angle in CM frame between the thrust axis of the 
rest of the event (ROE) and that of the B candidate is 
required to satisfy | cos(0tb.tb)| < 0.8 in order to re- 
duce background from e + e~ — > qq (q = u, d, s, c) con- 
tinuum events. The variable | cos($tb,tr)\ is strongly 
peaked near 1 for qq events, whereas BB events are 
more isotropic because the B mesons are produced 
close to the kinematic threshold. Additional separa- 
tion between signal and continuum events is obtained 
by combining several kinematic and topological vari- 
ables into a Fisher discriminant T, which we use in the 
maximum-likelihood fit described below. The variables 
|cos(07b,7h)|, \At\/a(Ai), |cos(0 B ,z)|, | cos(0 T s,z)|, and 
the output of a multivariate tagging algorithm [17] are 
used as inputs to T . The time interval At is calcu- 
lated from the measured separation distance Az between 
the decay vertices of B rec and the other B in the event 
(-Broe) along the beam axis (z). The vertex of T5 roc is 
reconstructed from the tracks that come from the sig- 
nal candidate; the vertex of t3r.oe is reconstructed from 
tracks in the ROE, with constraints from the beam spot 
location and the B ICC momentum. The uncertainty on 
the measured value of At is a (At). The variable 0b,z 
is the angle between the direction of B ICC and the z axis 
in the CM frame. This variable follows a sine squared 
distribution for BB events, whereas it is almost uniform 
for qq. The variable 6tb,z is the angle between the B 
thrust direction and the z axis in the laboratory frame. 

The decay modes studied are classified into three 
groups according to the final state particles: (i) B° — ► (fxfi, 
(ii) B+ -> 4>p+, and (iii) B° -» (j>p° , B° -> <f>f , and 
B° — ► /o/o- We find that 6% of events for the mode in 
group (ii) and 3% of events for the modes in group (iii) 
have more than one candidate that passes our selection 
criteria. For such events we retain the candidate with 
the smallest % 2 for the B rec vertex for use in the fits de- 
scribed below. The numbers of selected candidates are 
given in Table I. 

The dominant background for all modes comes from 
continuum events. The yield of this background compo- 
nent is determined from the fit to data. The dominant 
B backgrounds for group (i) are B° — ► (f)K*° and foK*°, 
which are estimated to contribute 1.4 and 0.6 events to 
the data, respectively. The B backgrounds for group (ii) 
are events from B decays to final states including charm 
and B + — ► (f>K* + . These are estimated to contribute 
107 and 5.5 events to the data. The B backgrounds for 
group (iii) are events from B decays to final states in- 
cluding charm, B° decays to 4>K*°, f K*°, 0/Y 2 *°(143O), 



and B + decays to (j)K + and <fiK* + estimated to con- 
tribute 249, 25.9, 9.1, 2.3, 4.7, and 1.8 events to the 
data. The branching fractions for the B backgrounds 
are taken from Ref. [18], except for B° — ► foK* , which 
has not yet been measured, and 4>p + where we use the 
results obtained here. The current upper limit on the 
B° — > foK*° branching fraction is 4.3xl0~ 6 and we as- 
sume a branching fraction of (2 ± 2)xl0~ 6 . 

We obtain yields for each mode from extended un- 
binned maximum likelihood (ML) fits with the input ob- 
servables Toes, AE, and cos^i^- In addition, for all 
modes except (fxfi, we include and T in the like- 

lihood, where mi^ is m^n or uikk for the <f>, p or fo 
candidates. A total of three fits are performed, one for 
each group of signal modes. We include event hypotheses 
for signal events and the aforementioned backgrounds in 
each of the fits. For each event i and hypothesis j, the 
likelihood function is 



i=l 



3=1 



where N is the number of input events, Nj is the number 
of hypotheses, rij is the number of events for hypothesis j 
and Vj (xj ) is the corresponding probability density func- 
tion (PDF) evaluated for the observables x, of the i th 
event. The correlations between input observables are 
small and are assumed to be negligible. Possible biases 
due to residual correlations are evaluated as described 
below. We compute the combined PDFs Vjfc) as the 
product of PDFs for each of the input observables. These 
combined PDFs are used in the fit to the data. 

For B decays to <p(j) and <j>p, the tttes distribution is 
parametrized with the sum of a Gaussian and a Gaus- 
sian with a low-side exponential component. The AE 
distribution is described by the sum of two Gaussian dis- 
tributions, and the cos #1,2 distributions are described by 
Eq. (1) multiplied by an acceptance function. The accep- 
tance function is a polynomial for all cos 01,2, with the 
exception of the p + helicity angle distribution for longi- 
tudinally polarized </>p + , which uses a polynomial multi- 
plied by the sigmoid function 1/(1 + cxp[a(cos 0i i2 + /?)]), 
where the parameters a and j3 are determined from MC 
simulated data. For the 4>p final states we use a Gaus- 
sian to describe the T distribution, and the sum of a 
relativistic Breit-Wigner (BW) with two Gaussians for 
7TT1.2- The continuum background toes distribution is de- 
scribed by an ARGUS function [19]. We parameterize the 
continuum AE distribution using a second-order polyno- 
mial and use polynomials to describe cos0i.2- Where 
appropriate, we parameterize the T distributions for the 
continuum background using a Gaussian, and we param- 
eterize the TO1.2 distributions using the sum of a BW 
and a polynomial. We use smoothed histograms of MC 
simulated data as the PDFs for all other signal and back- 
ground modes. We generate B° — > <fif assuming that the 
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TABLE I: Number of events N in the data sample, signal yield y*s (corrected for fit bias), fit bias, detection efficiency e, 
daughter branching fraction product (flft), significance a (including additive systematic uncertainties, taken to be zero if the 
fitted yield is negative), measured branching fraction where the first error is statistical, and the second systematic (see text), 
and the 90% C.L. upper limit on this branching fraction (including systematic uncertainties). For B decays to (f><j> and 4>p, two 
efficiencies are reported, one for longitudinally and one for transversely polarized events. The reported branching fractions for 
4>fo and /o/o are product branching fractions that are not corrected for the probability of fo decaying into 7r + 7r~ or K + K~ . 



Group 


N 


Mode 




Bias e(%) 


YlBi(%) 


a 


B( xlO -7 ) 


UL(xl0 -7 ) 


0) 


209 


<t>4> 


_1 r+3.7 
1 - J -2.9 


-0.4 ± 0.2 40.4 [28.7] 


24.3 ± 1.2 


0.0 


-0.4± -j ± 0.3 


<2.0 


(ii) 


3175 


4> P + 


oo c+11.3 


+2.3 ±1.1 5.7 [9.8] 


49.3 ±0.6 


2.2 


15± 7 6 ±9 


<30 


(iii) 


3949 


<t>P° 

<t>fo 
/o/o 


-i3.6±*:i 


+0.8 ±0.4 24.1 [26.5] 
-1.7 ±0.5 22.1 
-1.8 + 0.5 25.5 


49.3 ±0.6 


1.0 

0.0 
0.0 


0.9io.g ± 0.9 
0.2±° j ± 0.3 
— 1.4±g;| ± 1.5 


<3.3 
<3.8 
<2.3 



<f> is longitudinally polarised, and we use phase space dis- 
tributions for B° — » /o/o- Before fitting the data, we vali- 
date the fitting procedure using the methods described in 
Ref. [20]. We determine a bias correction on our ability 
to correctly determine the signal yield using ensembles 
of simulated experiments generated from samples of MC 
simulated data for the signal and exclusive backgrounds 
and from the PDFs for the other backgrounds. 

Our results are summarized in Table I where we show 
the measured yield, fit bias, efficiency, and the product of 
daughter branching fractions for each decay mode. We 
compute the branching fractions from the fitted signal 
event yields corrected for the fit bias, reconstruction ef- 
ficiency, daughter branching fractions, and the number 
of produced B mesons, assuming equal production rates 
of charged and neutral B pairs. As we do not know the 
value of Jl for the (jxp and <f)p modes, we fit the data 
for different physically allowed values of /z, in steps of 
0.1. We find no evidence for any of the signal modes and 
calculate 90% C.L. branching fraction upper limits xjjl 

such that $* UL L(y s , f L )dy s / j +o ° L(y s j L )dy s = 0.9, 

where L(y$, fa) is the likelihood as a function of signal 
yield 3^5 arid Jl multiplied by a uniform prior. We re- 
port the most conservative (largest) upper limits for each 
mode, for which //, = 0.5, 0.7, and 0.2 for groups (i), (ii), 
and (iii), respectively. The central values of the branch- 
ing fractions given in Table I correspond to these values 
of /l. Figure 1 shows the itjes distributions in subsam- 
ples of the data where \AE\ < 0.05 GeV for B+ -> </>p+, 
and \AE\ < 0.025 GeV for all other modes. 

We estimate the systematic uncertainty related to the 
parameterization of the PDF by varying each param- 
eter by its estimated uncertainty, and by substituting 
smoothed histograms by un-smoothed ones. The total 
contribution of all variations in signal yields, when added 
in quadrature, gives an error between 0.2 and 5.6 events, 
depending on the mode. We account for possible dif- 
ferences between data and MC events from studies of a 
control sample of B — > Dtt events, yielding an uncer- 
tainty of 0.1 to 12.2 events depending on the mode. The 




5.25 



5.26 



5.27 
n ES (GeV/c 2 ) 



5.28 



5.29 



FIG. 1: (color online) Signal-enhanced distributions of toes 
in data, with a projection of the fitted likelihood for (top) 
B° -» (j>cj>, (middle) B + -> 4>p + , and (bottom) B° -» 4>p°, 
B a — > </>/o, and B° — > /o/o- The solid line represents the 
total PDF, the dotted line represents signal, and the dashed 
line represents the sum of continuum and B backgrounds. 



uncertainty from fit bias is taken to be half the correc- 
tion listed in Table I. Incorporating the statistical un- 
certainty of the bias has a negligible effect. The uncer- 
tainty on _B-daughter branching fractions is in the range 
(1.2 to 4.9)% [18]. The modes in group (iii), <f>f , 
and /o/o have systematic uncertainties from the fo line- 
shape [21] of 0.2, 3.1, and 15.9 events, respectively. The 
mode B + — > </>p + has a fractional systematic uncertainty 
of 3.0% from the reconstruction efficiency of ir° mesons. 
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Other sources of systematic errors are track reconstruc- 
tion efficiency [(2.4 — 3.2)%], uncertainty on the number 
of B meson pairs (1.1%), particle identification efficiency 
(3.5%), and differences between data and MC efficiencies 
related to the cut on the vertex \ 2 (0.6%). 

Assuming isospin is conserved in / — > hh decays, 
where h = n,K, we correct for factors of B(fo — ► 
hh)/B(fo — ► h + h~) , to obtain the product branching 
fraction upper limits of B(B° — ► 0/ o ) x B(/o — > 7T7t) < 
5.7xl0" 7 , and £(B° - /„/„) x B(/ -» tttt) x S(/ - 
iiTA") < 6.9xl0~ 7 at 90% C.L. 

In summary we have performed searches for the decays 
B° — > 00, 0p°, 0/o, /o/o; an d -E^ - * 0P ± and place 
upper limits on these modes. The upper limit on B° — > 
00 reported here can be used to constrain possible NP 
enhancements suggested in Ref. [2]. 
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